Critical aging of a ferromagnetic system from a completely ordered state 
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We adapt the non-linear a model to study the nonequilibrium critical dynamics of O(n) sym- 
metric ferromagnetic system. Using the renormalization group analysis in d — 2 + e dimensions 
we investigate the pure relaxation of the system starting from a completely ordered state. We find 
that the average magnetization obeys the long-time scaling behavior almost immediately after the 
system starts to evolve while the correlation and response functions demonstrate scaling behavior 
which is typical for aging phenomena. The corresponding fluctuation-dissipation ratio is computed 
to first order in e and the relation between transverse and longitudinal fluctuations is discussed. 
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Aging phenomena have been found in a broad vari- 
ety of strongly disordered systems such as polymer and 
spin glasses [jj, electronic system of an Anderson in- 
sulator 0, array of flux lines pinned by disorder 
etc. Recently much attention has been attracted by ag- 
ing of pure (or weakly disordered) systems with slow 
relaxation dynamics governed, for example, by domain 
growth Q as in a ferromagnetic system below the crit- 
ical temperature T c or by critical slowing down as in 
a ferromagnetic system exactly at criticality 5J. One 
expects that the critical aging phenomena can be cast 
into different universality classes of nonequilibrium criti- 
cal dynamics. Most of studies consider the relaxation of 
a ferromagnetic system starting from a completely dis- 
ordered state after quenching it to the fixed tempera- 
ture T < T c . It was found that the response function 
R(t,s) and the correlation function C(t,s) depend non- 
trivially on the ratio x = t/s similar to that found in 
glassy systems. Here s and t are waiting and observation 
times respectively. The distance from equilibrium can 
be measured by the fluctuation-dissipation ratio (FDR) 
X(t,s) = TR(t, s) /d s C(t, s). It has been argued that 
for critical aging the limit X°° — lim^s—^ X(t. s) is a 
novel universal quantity of critical phenomena |5(. The 
FDR was computed for the d dimensional spherical model 
0, 0(n) symmetric ferromagnetic model @ and diluted 
spin models 0, H, @ ■ In all these systems X°° has values 
ranging between and 1/2. The mean field calculations 
show that the aging behavior is modified in the presence 
of long-range correlations in the initial disordered state 
[Tfl ITlj| . Much less known about the relaxation start- 
ing from an ordered state. The numerical simulations 
show that the correlation function demonstrates behav- 
ior which is typical for aging phenomena 0, ll2T | , while 
the magnetization obeys the long-time scaling behav- 
ior almost immediately after the system starts to evolve 
[l3l Il4| . This observation was used to develop new effec- 
tive numerical methods to determine the critical expo- 
nents, which are based on the short-time critical dynam- 
ics and do not require the time-consuming equilibration 
of the system |l5j . However, up to now there is no any 
theoretical explanation why the long-time scaling behav- 
ior emerges already in the macroscopically early initial 



stage of relaxation and there is no theoretical framework 
which allows one to take properly into account the critical 
fluctuations in this regime of aging. 

In this paper we study the nonequilibrium critical dy- 
namics of a ferromagnetic system starting from a com- 
pletely ordered state. The long-distance properties of the 
0(n) symmetric system below the transition point can be 
related to the non-linear a model defined by the reduced 
Hamiltonian 0] 



H = d d x 



i<v.)> 



h- s 



|s(x)| 2 = 1. (1) 



In principle, to describe properly the dynamics of a real 
isotropic magnet one has to consider the Larmor pre- 
cession of the conserved field s in the local magnetic 
field h, i. e. construct the low-temperature version of the 
model J [T^ . However in this work, which to our knowl- 
edge is the first analytical study of the critical aging of 
a finite range system from a completely ordered state, 
we are interested in the influence of the initial condition 
rather than in the complications caused by more realistic 
dynamic model. We restrict consideration to the pure 
relaxational dynamics (model A) with the non-conserved 
n-component order parameter s(x). It can be described 
by the Langevin equation 3 t s = — XSH/Ss + C, where A is 
an Onsager coefficient and £ is a random thermal noise 
[Til Il9ll20l| . Because of the constraint s 2 = 1 the random 
noise can act only in the tangential direction and thus is 
s - dependent. This dependence complicates the explicit 
expression for thermal noise distribution, however, one 
can simply formulate the dynamics in terms of a gen- 
erating functional |2(ij |. The corresponding weight reads 



P[s,s}= S{s 2 ~l)5(s-s)exp(-A), (2) 
where s is the response field and the action is given by 



A[s,s] = — I il''.nll ( As 2 



is ■ 



d t s + A— 

OS 



(3) 



Note that s is a real variable which is connected with the 
corresponding quantity s introduced in Refs. 



2 



by s = is. The first 5-function in Eq. (J2J is due to the 
constraint s 2 = 1 and the additional constraint s • s = 
0, imposed by the second 5-function, ensures that the 
thermal noise C acts only in the tangential direction [2(J . 

The effect of a macroscopic initial condition can be 
taken into account by averaging over the initial con- 
figurations Sq(x) = s(t — 0,x) with a weight P[sq] = 
exp(-H [s ]/T) |p. Taking Ho = - J d d x h • s (x) 
we specify an initial condition that corresponds to the 
equilibrium state of the system subjected to an exter- 
nal magnetic field ho, which is switched off at t = 0. 
There is a deep analogy between the considered model 
and the renormalization group (RG) description of the 
surface critical phenomena using expansion about the low 
critical dimension [2^. We assume that the magnetic 
fields ho and h act along the one direction which we will 
refer as to longitudinal one. We decompose the fields 
s = (<7, 7r) and s = (<r, #) into (n — f)-component trans- 
verse parts 7r and it and longitudinal parts a = \fl — tt 2 
and a = — tv ■ it /yl — tt 2 , where the constraints on the 
fields have been used. We now can write the weight func- 
tional exp(— At r [TT, 7r]) for the transverse components 
only with the action given by 

Ar[7r, tt] = ^ J d d x J dt {\-k 2 + iit(d t ~ AV 2 )tt 



+A 



(tt ■ tt) 2 

1-7T 2 



2 1 — 7T 



2\2 



(8 AV 2 W 2 A ( V?r 1 
(d t -XV )tt -^y^ 



+2\hVl-ir 2 ] } - j d d x^l~^ 



(4) 



We will adopt the dimensional regularization scheme [T^l 
in which the terms generated by the measure in the path 
integral over s and s vanish so that we have omitted these 
terms in Eq. from the beginning. Let us introduce the 
connected Green functions G"- :— ([n] k [ir] k ['Ko] l [iro] 1 ) , 
the low temperature expansions of which can be obtained 
with a loop expansion based on the action Q. The terms 
quadratic in tt and it give us the free response function 
and the free correlator: 



R°(t,s) 



C°Jt,s) = 



-i@{t - s)G q (t 
G g (\t-s\) 



h + q 



(5) 

2 ,G q (t + s ),(6) 



great advantage of the considered nonequilibrium model 
in comparison with the equilibrium counterpart is that 
the Green functions are not spoiled by ir singularities, 
so that we can probe the critical domain just taking the 
limit h — *■ without using the less convenient matching 
procedure [20l l23l | . However, the theory suffers of the uv 
divergences which can be converted into poles in e = d—2 
using dimensional regularization. Exploiting the ideas of 
Refs. ^(| and one can prove the renormalizability of 
the model, which means that all uv divergences can be 
absorbed into finite number of Z-factors according to 

7T = Z 1/2 7T, Tt = ZZ- 1/2 T\, 7T = (ZZ a ) 1/2 TV , 

i\ = ZZ-^Z^ttq, A = ZZ~ X \ f = yT e K d Z T T, 
h = n 2 Z T Z- 1/2 h, °h = [i 2 Z T {ZZoT lj2 ho. (7) 

Here circles denote the bare quantities, (2ir) d K,i is the 
surface area of a <i-dimcnsional unit sphere, and /i is an 
arbitrary momentum scale. The Z-factors except for Zq 
are the same as in equilibrium [Tsl Il9| and to order T 
given by Z= Z = l+(n-l)T/e and Z T = l + (n-2)T/e. 
The new factor Z = 1 + (n - 3)T/s + 0(T 2 ) serves to 
cancel the divergences arising from the free correlator © 
for t + s — > 0. The renormalized Green function reads 

Gf. =Z k l 2 {ZZ- l l 2 ) % {ZZ ) l ' 2 {Zl /2 ZZ- 1 / 2 ) l G l l % (8) 

and satisfies the RG equation 



^ + (3 T d T + (C - C)Ad A + phd h + { P + Co/2)h d ho 



+ (l + k-l-k)^ + (l + k)C + + l)y 



G^ = 0, (9) 



with /3 T = pdfj,T\ Q , Q = fid^ lnZ, ; | , p = f3 T /T + (/2-d, 
where | denotes the derivative at fixed bare parameters. 
Note that the equation similar to Eq. 10 holds for longi- 
tudinal and mixed Green functions. 

The critical behavior of the non-linear a model is con- 
trolled by the unstable ir fixed point (FP). For n > 2 
and d > 2 the beta function has a nontrivial zero 
T c = e/(n — 2) + 0(e 2 ), which is related to the bare 
critical temperature T c . The solution of Eq. 0, in con- 
junction with the simple dimensional analysis, yields the 
scaling behavior of the Green function at FP as 



where the notation G q (t) :— Te~ x ( h+q is introduced. 
The infinite number of higher order terms in the expan- 
sion of square roots in Eq. (@J will be treated as interac- 
tions. In each order in T we have to take into account 
only a fixed number of such terms. The completely or- 
dered initial state corresponds to the limit ho — > oo in 
which Eq. JHJ becomes a Dirichlet correlator. Although 
the action has infinite number of vertices proportional to 
ho, which are located at time surface t = 0, it is easy to 
show by direct inspection of Feynman diagrams that they 



do not contribute in the limit hn 



oo. The finite h n 



can be treated then as an additional perturbation. The 



G l l k (q,t;T,h,ho) = ^{Ty d ° M{T) l+k -" l -" k M{T) i+ " k 
xMo(T) l+i F l k [(q^ tC z ;hM£ d /T, h M /h) (10) 

where f(T) = pT x T x l e exp( J Q T dT'(l/(3 T (T') - 1/eT')) 
is the correlation length and da = d(k + k + I + I — 
1) — 2(1 + k) the canonical dimension of G l } t . The scal- 

1 v ' kk 

ing functions M, M and Mo are given by lnMi(T) = 
-\j T dT'Q(T')l(3 T {T'). Note that in statics M has 
a meaning of the spontaneous magnetization. For r = 
T - T c -> 0" we derive £ ~ \t\~ v and M - \rf, with 
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critical exponents v = —1/(3t(T c ) and f3 — is((T c )/2. 
The dynamic exponent reads z = 2 + C(T C ) — C(T C ). To 
one loop order we have v = 1/e, (3 = (n — l)/2(n — 2) 
and 2 = 2 + 0(e 2 ). 

We now focus on the scaling properties of two-times 
quantities which describe the critical system evolving 
from a completely ordered state. To that end we put 
in what follows h = 1 = and apply to Eq. IjlOjl a 
short-time expansion of the fields 7r and tt in terms of 
the initial fields -kq and 7To EH ■ As a result we ob- 
tain the scaling behavior of the response and correlation 
functions for s£~ z — > and t£~ z > 0: 



which is given by 



R q (t,s) = Ans-^+^/zit/s)* f R (\q z (t ~ s),t/s), (11) 
C q (t,s) = A c s^/ z (t/s) 5 f c (Xq z (t-s),t/s), (12) 



where 77 = 2/3/f — e is the Fisher exponent. The new 
dynamic critical exponent is defined as 9 — (2 — -q — 
z — (o(T c )/2)/ z, and should be distinguished from the 
initial slip exponent 9 |2l|. To lowest order in e it is 
given by 9 — —e(n — l)/4(n — 2) + 0(e 2 ) and coincides 
with —(3/vz. Although at this point we do not see any 
symmetry which ensures this identity, we argue that it 
holds to all orders in e and below give some arguments 
supporting this conjecture. The functions f R and fc are 
regular functions of both arguments and finite for s — > 0. 
We have explicitly computed the renormalized response 
and correlation functions to one-loop order: 

R q (t,s) = K- 1 n[t, S }R° q (t,s), (13) 
C q (t,s) = K^[t, S ]C° q {t,s) 

+ (n-2)T 2 e- Xq2( - t+s) F{2Xq 2 s)/2K d q 2 , (14) 

where we have introduced notations Q[t, s] = 1 — ^(7 + 
ln2A/x 2 s+2=ilnf) and F(y) = 7 + ln \y\ -Ei(y). Here 7 
is the Euler constant and Ei(cc) the exponential integral, 
so that F(0) = and F'(0) = -1. Substituting the FP 
value T c in Eqs. I|13|) and l|14|l we find that they are consis- 
tent with scaling laws i|ll|) and (|12|) . The corresponding 
non-universal amplitudes to one loop order are given by 
A R = K-\-iT c ){l - 7T c /2)(2Am 2 )- t ^ 2 , A c = i\A R , 
and the universal (apart from the normalization) scaling 
functions to the same order read 



fn{u,v) = e u , 
fc(u,v) = (v - l){e" 



-u(v+l)/(v-l) 



x [l-eF(2u/(v-l))/2]}/u. 



(15) 



(16) 



Let us introduce the FDR for a particular mode q as 6] 
X- 1 = d s C q (t,s)/i\R q (t,s). Using Eqs. (O and dJ| 
we obtain to one loop order that X" 1 = fx(2Xq 2 s) with 

fx(u) = l + e- u + £ [(n-3)(l-e-")/(n-2)u 

- e- u (F(u)-2F'(u))}/2. (17) 

According to arguments of Rcf. 6] the local FDR X°° is 
identical to the global FDR = lim^s^oo X q= o(t, s), 



XZo = tt(1 



In- 1 



4)i-2 



0(e 2 ). 



(18) 



Note that X^ > ^ while for most systems, except for 
some urn models [23, the aging from a disordered state is 
characterized by the FDR X°° S [0, |]. Thus the system 
aging from an ordered state is more close to equilibrium 
(X = 1) than when it relaxes from a disordered state. 
We expect that this applies not only to the system under 
consideration but also to all systems with pure relaxation 
dynamics including the Ising model. 

We are now in a position to discuss the relaxation of the 
magnetization. The asymptotic long-time behavior of the 
magnetization can be deduced directly from Eq. I|1U|) for 
Giq and reads M(t) cx t~®l vz . To describe the behavior 
of the magnetization during the macroscopically early ini- 
tial stage of relaxation we derive the equation of motion 
for the magnetization. The latter can be written in the 
form of a Ward identity J °° ds T q= o(t, s)M(s) = 0, which 
follows from the invariance of the generating functional 
under simultaneous rotation of fields s and s [l9t |20| . 
Here T q (t, s) is a vertex function which can be expressed 
as a sum of one-particle irreducible diagrams with am- 
putated external lines 7r and 7r. To one- loop order we 
have T q {t, s) = iT- 1 6(t - s)V s (q) with the renormalized 
operator V t — ZZ^}T> t given for q = by 

V t (0) = d t + ^ [7 + ln2A^ 2 t] d t +^-j+0(T 2 ). (19) 

The equation of motion can be written as T> t (Q)M(t) = 0, 
t > and T ~T C . Solving this equation to first order in e 
we have to neglect the second time derivative in Eq. H19|) 
since it carries a factor T and hence in the equation of 
motion only gives rise to a correction of order e 2 . Having 
done that, we get M(t)_= Mt^t"- 1 )/^"- 2 ), where the 
non-universal constant M is determined by motion of the 
system during the microscopically early initial stage of 
relaxation. Indeed, the considered field theory describes 
the behavior of a ferromagnet only on scales larger than 
the lattice constant and on times larger than the corre- 
sponding microscopic time T m ; c . Therefore, the derived 
power-law behavior, which is usually called the non- linear 
critical relaxation, does exist only for r m i C < t < r mac . 
Here T mac is the macroscopic time at which there is a 
crossover to linear relaxation M{t) cx exp(— i/r mac ). It is 
determined either by size of the system r mac cx L z or by 
temperature r mac cx \t\~ zv . 

Let us show how this picture can be generalized to 
higher orders. The vertex function is related to the cor- 
responding Green function by the Dyson equation [l9l | 
L°° dt'Y q {t,t')R q {t' ', s)(s) = S(t — s), which is written in 
the time domain and can be easily check to one loop by 
means of Eqs. (1K'>[1 and l|19|) ■ Using Eq. I|ll|) we obtain 

r q=0 (t,s)=A^- 2 +^ z {[S(t-s) 

- (t/s) S f'(t/s)/s\ (d s -9/ S ) + (t/ s fyj(t/ s )/s 2 }, 
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where f(v) = fn(0,v) and -0 is the solution of Voltera 
equation £ du[f'(u)f'(v/u)v/u — f(v/u)ip(u)] = such 
that = 0. This implies that the solution of the 

equation of motion can be written as M (t) = Mt e fM(t) 
with fmit) satisfying the condition /m(oo) = 1. Taking 
into account the known asymptotic behavior for t — > oo 
we conclude that the identity 9 — —f3/vz holds to all 
orders. To find the exact form of /m we have to solve 
the corresponding integral equation. Recently, the local 
scale invariance (LSI) was used to predict the exact scal- 
ing form of the response function for a system evolving 
from a completely disordered state j^. We expect that 
the transverse fluctuations of the considered model share 
similar scaling properties. In analogy with Ref. |26j| we 
may expect that f(v) = 1 and consequently ip(v) = and 
/m(£) = 1 to all orders. However, the analysis of two loop 
graphs suggests that in accordance with Ref. there 
should be small corrections to LSI: f(v) = 1 + 0(e 2 ). 
They come from the non-local in time contributions to 
the vertex function T q (t, s) and give rise to very small 
corrections in the magnetization /m(^) = 1 + 0(e 4 ) that 
are likely difficult to observe in simulations. 

We now discuss the relation between trans- 
verse and longitudinal fluctuations. Let us 
define the longitudinal response function as 
K q (t,s) = Jd d x e lqx (cr(x, t)a(0, s)) c and analo- 
gously the longitudinal correlation function C q (t,s). 
Expressing a and a in terms of 7r and ir we obtain 
their low temperature expansions. To lowest order the 
renormalizcd longitudinal response function is given by 



K q (t,s) 



2K d 



-T 



In- 



t - s 



+ F(Xq 2 (t~s) 2 /2t) 



F (Xq 2 (t - s)/2)] R% s) + 0(T 2 s, T 3 ). (20) 



The corresponding expression for the longitudinal corre- 
lator is too cumbersome so that for the sake of conciseness 
we write down only its time derivative at q = 



d s C q=a {t,s) 



2K d 



VAln t + S 



t 



0(T 2 e,T 3 



(21) 



which is needed to compute the longitudinal FDR X q = 
i\lZ q (t,s)/d s C q (t,s). For the latter we obtain X£1 Q = 
\ + 0(e). The computation to order e requires two-loop 
calculations and is left for future investigations. However, 
we expect that the scaling ansatzs and (I12f> are valid 
also for the longitudinal functions. 

After we submitted our preprint to arxiv.org there 
appeared Refs. an d HH which also study the crit- 



ical aging from a magnetized state but for infinite range 
models such as the spherical model, long range ferro- 
magnetic model and Ising model in the limit of large 
dimension. These works confirm that the critical aging 
from an ordered state is characterized by A°° > 1/2. 
The d — 2 expansions of the response and correlation 
functions obtained in Ref. [2^| for the spherical model 
are in agreement (up to prefactors) with our results for 
the longitudinal functions l|20|) and i|21|) . The difference 
between the spherical limit (n — ► oo) of the transverse 
FDR (O and the FDR of the spherical model [H sug- 
gests that X°° ^ X°° for finite n and thus one cannot 
introduce the observable-independent effective tempera- 
ture T c ff = T/X°° 30]. Rewriting the response func- 
tion of the spherical model in the form i|ll|) we obtain 
/^P(0,u) = 1 - (1 - l/vy d ~ 2 ^ 2 . Thus in contrast to the 
transverse response of the 0(n) system with the scal- 
ing function fn(0,v) being finite in the limit v — * oo, 
the scaling function of the spherical model decays as 
/|j P (0, v) oc v~ x . Supposing the same asymptotics for the 
longitudinal function fn(0,v) of the O(n) system we ar- 
rive at 7^=o (*, s) oc s-O-^/^t/sf _1 for t > s. This, 
of course, has to be checked independently. 

In summary, we have studied the nonequilibrium crit- 
ical behavior of a ferromagnetic systems which evolves 
from a completely ordered state. We have shown that 
the universal long-time scaling behavior of magnetiza- 
tion emerges in the macroscopically early initial stage of 
relaxation and that can be considered as a consequence 
of the LSI. The two-times functions exhibit aging behav- 
ior which in contrast to the aging from a disordered state 
is characterized by the FDR larger than 1/2. Numerical 
simulations of Ref. 0] show that the generalized scaling 
behavior exists in the short-time regime of the critical 
relaxation even for an initial state with arbitrary magne- 
tization. However, instead of the critical exponent 9 there 
exists a whole universal characteristic function which up 
to now can be estimated only numerically. We expect 
that our computation being extended to finite ho can de- 
scribe this scaling behavior and provide a way to calcu- 
late the corresponding characteristic function. We hope 
that the considered model can also be useful to study the 
influence of Goldstone modes on phase ordering kinetics 
at T < T c . 
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